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Introduction 

The  use  of  lasers  in  Air  Force  (AF)  operations  is  increasing  rapidly  (Scott,  1999), 
and  the  deployment  of  low-power  lasers  as  anti-personnel  weapons  is  a  growing  threat  to 
AF  aircrews  (Mellerio,  Marshall,  Tengroth,  Anderberg,  &  Wolbarsht,  1991). 

Consequently,  the  possibility  of  a  loss  in  vision  resulting  from  an  accidental  or  hostile 
laser  exposure  continues  to  increase  year  after  year.  The  Air  Force  Research  Laboratory 
Optical  Radiation  Branch  (AFRL/HEDO)  has  been  studying  this  problem  for  over  two 
decades.  This  effort  has  led  to  the  creation  of  mathematical  models  that  describe  how 
basic  visual  processes  are  affected  by  exposure  to  optical  sources,  such  as  lasers.  Such 
models  are  important  because  they  help  analysts  understand  how  a  functional  visual 
deficit  impacts  aircrew  performance  and  ultimately  mission  success. 

Exposure  to  a  visible  laser  beam  can  cause  temporary  visual  impairment,  such  as 
glare  and  flashblindness,  without  producing  physical  eye  damage.  AFRL/HEDO  has 
developed  models  of  these  effects.  The  Glare  Model  is  principally  based  on  the  work  of 
Vos  and  others  at  the  TNO  Human  Factors  Research  Institute  (Vos,  1984;  Vos  &  van  den 
Berg,  1997).  The  Flashblindness  Model  was  originally  developed  by  Menendez  and 
Smith  (1990)  at  AFRL/HEDO.  Though  related,  the  models  are  distinguished  by  the 
timing  of  the  optical  radiation  source  -  the  Glare  Model  estimates  the  visual  impairment 
during  the  laser  exposure,  while  the  Flashblindness  Model  predicts  residual  effects  after 
the  exposure. 

This  report  reviews  the  state  of  glare  and  flashblindness  modeling  at 
AFRL/HEDO.  The  scientific  basis  of  the  algorithms  is  presented  and  the  appropriate 
research  literature  is  reviewed.  The  strengths  and  weaknesses  of  the  approaches  are 
compared  and  evaluated  in  terms  of  their  applicability  to  AF  optical  radiation  concerns.  In 
addition,  issues  that  relate  to  the  assumptions,  validity,  and  scope  of  the  models  are 
addressed.  These  issues  include  short  pulse  exposures,  multiple  pulse  exposures,  and  the 
impact  of  exposure  intensity  and  duration  on  light  adaptation. 

The  Glare  Model  is  presented  first.  The  major  components  of  the  model  are 
presented,  including  optical  radiation  source  parameters,  transmission  of  light  through  the 
eye,  the  principle  of  the  equivalent  background,  the  concept  of  veiling  glare,  and  target 
considerations.  The  report  continues  with  a  discussion  of  the  Flashblindness  Model.  The 
Flashblindness  Model  includes  all  of  the  components  in  the  Glare  Model  with  the  addition 
of  a  flashblindness  recovery  function  that  describes  the  process  of  dark  adaptation.  Three 
candidate  algorithms  for  the  flashblindness  recovery  function  are  discussed  and  compared. 
The  current  model  (Menendez  &  Smith,  1990)  is  compared  to  two  alternative  models  -  an 
updated  one  by  Smith  (1996)  and  one  adapted  from  Hahn  and  Geisler  (1995).  The  model 
that  best  represented  the  experimental  flashblindness  recovery  data  was  incorporated  in 
LTMC.  Finally,  recommendations  to  further  develop  and  improve  the  Glare  and 
Flashblindness  Models  are  made. 
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Laser  Threat  Modeling  Component  Application 

The  Flashblindness  and  Glare  Models  are  parts  of  the  Laser  Threat  Modeling 
Component  (LTMC)  application.  LTMC  is  the  flagship  laser  hazard  bioeffects  model  of 
AFRL/HEDO.  LTMC  is  a  Microsoft  Windows  NT™-based  integrated  modeling  and 
simulation  environment  that  helps  analysts  evaluate  the  vulnerability  of  aircrew  to  optical 
radiation  threats  in  operational  and  training  environments.  LTMC  is  composed  of  a  set  of 
modular  software  components  that  provides  hazard  evaluations,  classification  and  range 
calculations  for  optical  radiation  bioeffects  including  eye  damage,  flashblindness,  and 
glare.  The  Flashblindness  and  Glare  Models  provide  estimates  of  reductions  in  visual 
sensitivity  as  a  function  of  distance  from  the  optical  source.  This  information  allows 
mission  planners  to  assess  aircrew  vulnerability  during  a  mission  profile. 

LTMC  Glare  Model 

The  Glare  Model  determines  how  the  visibility  of  a  target  is  degraded  by  glare. 

For  this  analysis  the  model  uses  quantitative  measures  of  the  source  exposure,  the  target 
and  background  scene  brightness,  and  the  visibility  threshold  of  the  target  of  interest.  The 
model  first  converts  the  source  exposure  from  radiometric  units  into  photometric  units. 

An  intraocular  scatter  function  describes  how  the  light  source  is  distributed  on  the  retina 
after  passing  through  the  ocular  media.  The  scatter  function  transforms  the  source  into  a 
glare  field  of  a  given  luminance  distribution.  The  glare  field  produces  a  relative  scotoma, 
that  is,  an  area  of  reduced  visual  sensitivity  in  the  field  of  view.  The  scotoma  is  formed  by 
superimposing  the  glare  field  onto  the  target  and  background  scene,  which  reduces  scene 
contrast.  The  model  compares  the  reduced  target  contrast  to  the  target  visibility  threshold 
and  determines  if  the  target  has  been  obscured.  A  detailed  description  of  the  Glare  Model 
follows  beginning  with  the  radiation  source. 

Photometric  quantities 

Although  laser  sources  are  usually  described  in  terms  of  radiometric  quantities, 
these  give  no  indication  of  the  effectiveness  of  the  source  as  a  stimulus  for  vision.  The 
photometric  system  was  developed  to  describe  optical  radiation  in  terms  of  its  ability  to 
elicit  a  response  from  the  visual  system.  The  basic  photometric  unit  of  luminous  power, 
or  luminous  flux,  is  the  lumen  (lm).  The  lumen  is  the  photometric  equivalent  of  the  watt 
(W).  Under  photopic  conditions,  one  watt  of  radiant  flux  at  555  nm  is,  by  definition, 
equivalent  to  685  lm  of  luminous  flux.  This  wavelength  is  at  the  peak  of  the  spectral 
sensitivity  of  the  eye  under  photopic  viewing  conditions.  The  relative  sensitivity  of  the 
eye  is  defined  by  the  Commission  Internationale  de  V Eclair age  ( CIE )  (1926)  relative 
photopic1  spectral  luminosity  function,  Vx ,  for  a  standard  observer  (Figure  1).  This 
function  is  given  a  value  of  unity  at  555  nm,  its  maximum. 

For  a  laser  with  a  radiant  flux  of  P  (in  W)  at  555  nm,  the  luminous  flux  is 
685  x  P  lm.  To  generalize  to  any  laser  wavelength,  X,  the  luminous  flux,  is  calculated 
as: 


1  The  photopic  function  is  used  because  photopic,  or  cone  vision,  is  used  for  essentially  all  aircrew 
visual  functions. 
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$  =  685xPxVA9  (1) 

where  P  is  the  radiant  flux,  and  Vx  is  the  value  of  the  relative  spectral  luminosity  function 
at  this  wavelength. 
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Figure  1,  The  CIE  photopic  spectral  luminosity  curve  for  the  standard  observer 
(1926).  The  curve  shows  peak  sensitivity  at  555  nm  and  declines  toward  zero  at  400 
and  750  nm. 


Optical  source  considerations 

The  Glare  Model  analyzes  glare  effects  for  continuous  wave  (CW),  single 
wavelength,  point  source  exposures,  based  on  the  source  irradiance  measured  at  the 
cornea.  The  photometric  equivalent  of  irradiance  is  illuminance,  Ev  (in  limin'2  m  lux), 
which  is  related  to  irradiance,  Ee  (in  W-cm  "2),  by: 

Ev  =  685  xEexVxx  104 .  (2) 

In  addition,  since  the  human  visual  system  perceives  a  flickering  light  source  with 
a  frequency  above  60  Hz  as  continuous  (Brindley,  1960),  the  effects  of  pulsed  laser 
systems  with  pulse  repetition  frequencies  greater  than  60  Hz  can  also  be  analyzed.  For 
these  sources,  the  irradiance  is  given  by  the  time-averaged  irradiance. 

A  schematic  view  of  the  eye  is  shown  in  Figure  2.  Light  rays  from  a  distant 
optical  source  are  focused  on  the  retina  and  are  concentrated  approximately  50,000  times 
(Sliney  &  Wolbarsht,  1976).  This  is  what  makes  a  visible  laser  so  potentially  damaging  to 
ocular  tissue  compared  to  skin  tissue.  Most  lasers,  when  viewed  from  a  distance,  are 
considered  point  sources,  because  they  appear  to  have  little  or  no  spatial  extent.  The 
American  National  Standards  Institute  (ANSI)  provides  guidelines  for  determining 
whether  a  laser  may  be  considered  a  point  or  an  extended  source  (American  National 
Standards  Institute,  2000).  The  standard  defines  a  laser  that  subtends  an  angle  less  than 
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some  minimum  angle  at  the  eye  as  a  point  source.  The  criterion  minimum  angle  depends 
on  exposure  duration,  but  when  the  source  produces  an  image  on  the  retina  of  20-30  pm  in 
diameter,  and  is  nearly  diffraction-limited,  then  it  is  generally  considered  a  point  source. 


Figure  2.  A  schematic  diagram  of  the  human  eye  showing  how  light  from  a  distant 
source  forms  an  image  on  the  retina.  If  the  source  subtends  an  angle  less  than  the 
critical  angle,  at,  then  it  is  considered  a  point  source. 


Retinal  distribution  of  a  point  source 

Although  the  light  from  a  point  source  should  be  imaged  on  the  retina  as  a  point,  in 
reality,  the  focused  rays  do  not  all  converge  on  one  spot.  Optical  imperfections  such  as 
spherical  and  chromatic  aberrations,  diffraction  effects,  and  refractive  errors,  cause  the 
image  of  a  point  of  light  to  be  blurred  on  the  retina  (Campbell  &  Green,  1966;  Campbell 
&  Gubisch,  1966;  Vos,  1984;  Westheimer,  1986).  Retinal  image  quality  is  further 
degraded  by  intraocular  stray  light  from  internal  ocular  reflections. 

Estimates  of  the  extent  of  laser  glare,  therefore,  must  begin  with  a  calculation  of 
how  the  laser  light  is  distributed  on  the  retina.  The  precise  shape  of  the  retinal  distribution 
has  been  a  topic  of  great  interest  in  applied  vision  work  related  to  veiling  glare  and  there 
has  been  a  concerted  effort  to  quantify  it.  A  number  of  studies  have  been  compiled  and 
analyzed  by  Vos  and  Berg  (1997)  who  have  produced  an  analytical  expression  for  the 
point  spread  function  (PSF)  that  closely  fits  most  available  data,  distinguishes  the  different 
components  of  entoptic  scatter,  and  explicitly  depends  on  age  and  complexion.  The 
expression  breaks  the  sources  of  intraocular  scatter  into  three  components:  an  anterior 
component  from  the  cornea  and  lens;  a  wall  component  from  the  iris  and  sclera;  and  an 
epithelium  component  for  scatter  from  the  retinal  pigment  epithelium  (RPE).  The  PSF 
describes  a  radially  symmetric  distribution  of  light  on  the  retina  as  a  function  of  angular 
distance,  0(in  degrees),  from  the  center  of  the  source: 

PSF  (9)  =  Anterior(0)  +  Wall(0)  +  Epithelium (#)  (3) 

The  expression  for  scatter  by  the  anterior  segment  contains  an  age  factor  (Ijspeert, 
de  Waard,  van  den  Berg,  &  de  Jong,  1990;  Ijspeert,  van  den  Berg,  &  Spekreijse,  1993;  van 
Meeteren,  1974)  and  is  given  by: 
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Anterior(8)  =  1.03x 


where  A  is  the  person’s  age  in  years  and  0is  the  glare  angle  in  degrees. 

The  amount  of  light  passing  through  the  iris  and  sclera  is  given  by; 


Wall(9)=  3.5  xlO"5 


(105 -0)  ^ 

eos(O.92(0-5)), 


(5) 


The  magnitude  of  the  epithelial  is  component  depends  on  the  degree  of 
pigmentation,  with  a  heavily  pigmented  RPE  scattering  less  light  than  a  lightly  pigmented 
one,  and  is  given  by: 


Epithelium^  =  0.17  5P 
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where  P  is  the  pigmentation  of  the  eye,  which  varies  from  0  in  heavily  pigmented  non- 
Caucasians,  to  1.2  for  very  light-eyed  Caucasians  (van  den  Berg,  1995).  The  variation  in 
the  PSF  due  to  age  (in  the  typical  range  for  aircrew)  and  pigmentation  is  small,  so  a 
general  function  is  used  in  die  Glare  Model  (Figure  3),  which  is  based  on  a  30  year  old, 
brown-eyed  Caucasian  (P=0.75). 

Another  algorithm  has  been  developed  that  estimates  the  retinal  light  distribution 
from  an  extended  light  source  (Toet,  Ijspeert,  Vos,  &  Walraven,  1995).  This  algorithm 
convolves  the  PSF  with  an  extended  light  source  to  predict  a  two-dimensional  retinal  light 
distribution.  When  added  to  LTMC,  this  algorithm  will  be  able  to  predict  the  effects  of 
extended  light  sources,  such  as  exposures  passing  through  optical  media  (e.g.  aircraft 
canopies). 

It  is  known  that  intraocular  scatter  increases  when  light  enters  the  eye  obliquely, 
that  is,  not  parallel  to  the  visual  axis.  Oblique  scatter  increases  appreciably  when  the 
angle  of  incidence  exceeds  about  20°  (Jennings  &  Charman,  1981).  However,  because  the 
Glare  Model  assumes  the  point  source  is  always  parallel  to  the  visual  axis,  no  corrections 
are  necessary  for  obliquely  incident  light. 

Since  the  Glare  Model  is  a  static  model,  it  assumes  that  the  glare  source  and  target 
are  fixed  in  space  and  time.  It  does  not  account  for  eye  movements,  which  means  that 
glare  effects  are  fixed  with  respect  to  retinal  location. 


The  equivalent  background  luminance 

The  effort  to  measure  the  effect  of  a  glare  source  on  the  retina  has  a  long  history 
(Vos,  1963, 1984).  In  an  early  paper  Cobb  (191 1)  introduced  the  concept  of  the 
equivalent  background  luminance  (EBL)  to  extrapolate  the  effect  of  a  glare  source  to  a 
wide  variety  of  visual  tasks.  He  showed  that  vision  in  the  presence  of  a  glare  source  was 


5 


impaired  as  if  a  veil  of  light  were  cast  over  the  targets  in  the  field  of  view.  The  visual 
effect  of  the  source  could  be  described  by  equating  the  glare  intensity  at  an  angular 
distance  on  the  retina  to  a  uniform  distribution  of  light  in  visual  space.  The  two  were  said 
to  be  equal  when  both  sources  had  the  same  masking  effect  on  a  target  of  a  specific 


Figure  3.  The  point  spread  function  describing  intraocular  light  scatter  in  the 
human  eye.  The  PSF  is  radially  symmetric  about  the  visual  axis  (0°)  and  is  shown 
here  for  a  30-year  old  brown-eyed  Caucasian. 

luminance.  Thus,  the  intensity  of  a  glare  source  at  any  point  along  its  profile  could  be 
expressed  as  a  quantifiable  light  veil  external  to  the  eye.  This  comparison  was  performed 
psychophysically,  where  an  observer  looked  at  a  point  a  certain  angle  from  the  center  of 
the  glare  source  and  adjusted  the  brightness  of  an  extended  light  source  that  had  a 
comparable  masking  effect  on  a  given  target.  This  process  was  repeated  over  a  range  of 
angles  to  provide  an  EBL  distribution  of  the  glare  source2.  The  results  were  usually 
expressed  as  a  PSF  by  plotting  the  equivalent  luminance  of  the  glare,  Leg,  relative  to  the 
incident  illumination  from  the  glare  source,  Eg?. 

PSF(&)  =  Sl  (7) 

where  is  given  in  cd-m'2  and  the  illuminance  of  Eg!  is  in  lm  m‘2. 


2  The  fact  that  the  ratio  of  the  glare  source  intensity  to  the  equivalent  background  luminance  was 
constant  over  a  wide  range  of  intensities  supports  the  argument  that  the  source  of  the  glare  was  scatter  from 
optical  surfaces,  and  not  neural  in  nature  (Vos,  1984). 

3  In  photometric  units,  this  is  expressed  in  units  of  sr'1  (derived  from  zA  rci2  per  lux,  which  is 
equivalent  to  lm-sr^m"2  per  lm-m"*). 
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By  simple  rearrangement  of  Equation  (7),  the  retinal  distribution  of  a  glare  source 
can  be  expressed  as  an  EBL  profile  by  multiplying  the  glare  illumination  at  the  cornea  by 
the  PSF: 


Leq  =  Egl  x  PSF(0) ,  (8) 

The  advantage  of  stating  the  glare  source  in  terms  of  an  EBL  is  that  the  EBL  can  be 
expressed  in  the  same  brightness  units  (e.g,,  cd-m'2)  as  targets  and  backgrounds.  In  this 
way  glare  can  simply  be  added  to  the  target  and  background  to  determine  its  effect  on 
target  visibility. 

The  Glare  Model  limits  the  glare  intensity  as  well.  An  upper  limit  is  based  on  the 
premise  that  if  all  of  the  photopigment  is  bleached  by  a  given  exposure,  then  any 
additional  light  has  no  effect  on  the  light  adapted  state  of  the  eye.  There  is  general 
agreement  that  a  7.6  log  td4  exposure  bleaches  100%  of  the  photopigment  (Rushton  & 
Henry,  1968).  Any  additional  photons  are  not  absorbed  and,  therefore,  do  not  result  in  a 
change  in  sensitivity  or  in  light  adaptation.  Thus,  7.6  log  td  was  set  as  the  limit  for 
modeling  glare  intensity5.  The  next  section  discusses  the  way  target  visibility  is  reduced 
by  glare. 


Computing  glare  effects  on  target  visibility 

Once  the  light  distribution  of  the  glare  source  on  the  retina  is  computed  and 
expressed  in  terms  of  an  EBL,  then  the  effect  on  target  contrast  can  be  estimated. 

The  contrast,  C,  of  a  target  in  the  visual  scene  is  computed  as: 


C  = 


(9) 


where  T  is  the  target  luminance  and  B  is  die  background  luminance,  both  of  which  are 
expressed  in  cd-m"2.  The  absolute  value  of  the  difference  between  the  target  and 
background  luminance  is  taken  because  the  difference  in  the  visibility  of  positive  contrast 
and  negative  contrast  targets  is  negligible. 


The  superimposition  of  the  veiling  glare  source  is  expressed  mathematically  by 
adding  the  EBL  to  the  target  luminance  and  background  luminance  quantities.  The 
apparent  target  contrast,  C\  in  the  presence  of  the  glare  source  is  given  as: 


| (T  +  EBL)- (B  +  EBL) |  _  | T-B 
B  +  EBL  ~  bTebE’ 


(10) 


The  veiling  luminance  from  the  glare  source  does  not,  therefore,  affect  the 
difference  between  the  target  luminance  and  the  background  luminance,  but  provides  an 
elevated  background  against  which  this  luminance  difference  must  be  viewed  and,  hence, 
reduces  target  contrast.  Since  the  glare  decreases  with  angular  distance  from  the  center  of 


pupil 


4  A  troland  is  defined  as  the  retinal  stimulation  provided  by  a  source  of  1  cd-m"2  viewed  through  a 
of  area  of  1  mm2. 

5  The  Glare  Model  assumes  a  one-second  exposure. 
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the  image  of  the  source,  that  is,  EBL  is  a  function  of  6  (Equation  8),  the  apparent  target 
contrast  must  be  computed  at  each  retinal  location  across  the  glare  profile.  The  Glare 
Model  assumes  that  the  target  and  the  veiling  glare  are  radially  symmetric  so  the  apparent 
target  contrast  is  only  computed  along  one  meridian. 

To  determine  whether  the  glare  source  obscured  the  target,  the  apparent  target 
contrast  is  compared  to  the  target  contrast  threshold.  The  procedure  must  be  carried  out 
for  each  point  along  the  retinal  meridian.  If  the  apparent  target  contrast  is  below  the 
contrast  threshold,  then  the  target  is  not  visible  at  that  point.  The  target  is  considered 
obscured  when  the  target  is  below  threshold  across  its  entire  extent. 


Target  considerations 

The  Glare  Model  assumes  that  the  observer  is  looking  directly  at  the  target  and  that 
the  optical  source  is  close  enough  to  the  target  so  that  the  images  are  effectively  co¬ 
incident.  Thus,  both  the  target  and  the  peak  of  the  glare  source  are  centered  on  the  fovea. 
This  arrangement  yields  the  maximum  likelihood  of  obscuring  the  target  and  represents  a 
worst-case  scenario  for  target  detection. 


Targets  are  expected  to  be  small  in  angular  extent,  because  they  are  usually  viewed 
from  long  range.  Therefore,  target  detection  algorithms  that  pertain  to  small,  spot  targets 
are  appropriate  for  use  in  the  model.  The  default  target  detection  algorithm  is  from  the 
visual  detection  model,  VIDEM,  a  model  created  for  predicting  aircraft  detectability 
(Akerman  &  Kinzly,  1979).  The  algorithm  estimates  the  contrast  necessary  for  the  target 
to  be  just  visible  (50%  probability  detection  criterion)  against  an  unstructured,  uniform 
background.  The  algorithm  is  based  on  a  set  of  laboratory  contrast  threshold  tests  (Lamar, 
Hecht,  Schlaer,  &  Hendley,  1947, 1948;  Sloan,  1961),  which  were  then  validated  by 
aircraft  detection  field  tests.  Using  these  data  Akerman  and  Kinzly  (1979)  found  that  the 
contrast  threshold,  Ct,  was  best  described  by  the  following  function: 


C, 


0.0352  x  0024 


0.584  xG16 


(11) 


where  0is  the  distance  of  the  target  from  the  fovea  center  (in  degrees)  and  a  is  the  target 
size  (in  minutes-of-arc). 

Figure  4  shows  a  family  of  VIDEM  contrast  thresholds  for  various  target  sizes  and 
angular  distances  from  the  fovea.  The  VIDEM  algorithm  was  developed  to  predict  the 
contrast  threshold  for  detection  of  aircraft  against  daytime  sky  background  using  aircraft 
size  and  retinal  location  as  variables  and  assumes  that  the  target  is  small  and  circular  in 
shape.  If  the  target  is  not  circular,  then  the  target  must  be  expressed  in  terms  of  an 
equivalent  circular  area. 

One  limitation  of  VIDEM  is  that  the  peripheral  threshold  data  are  all  based  on  a 
single  level  of  background  adaptation  (10  cd-m-2).  This  level  would  be  appropriate  for 
detecting  aircraft  just  after  dawn  and  just  before  dusk.  Indeed,  increment-detection 
measures  indicate  that  Weber’s  law  of  constant  thresholds  begins  around  10  cd  m"2  for 
small  spot  targets  (Blick,  Beer,  Kosnik,  Troxel,  Toet,  Walraven,  &  Mitchell,  2001 .  This 
implies  that  contrast  thresholds  are  the  same  at  higher  adaptation  levels  and  no  adjustment 
for  higher  background  luminance  levels  is  necessary.  However,  contrast  thresholds  are 
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higher  at  luminance  levels  below  10  cd-m'2  (Smith,  1996),  and  this  argues  for  an 
expansion  of  the  VIDEM  database  if  detection  thresholds  for  lower  adaptation  levels  need 
to  be  estimated. 


Figure  4,  The  VIDEM  model  predicts  the  detectability  of  aircraft  against  a  daytime 
sky  background  (Akerman  &  Kinzly,  1979).  Contrast  thresholds  are  based  on 
target  diameter  in  minutes  of  arc  and  retinal  location  in  degrees.  The  model 
assumes  radial  symmetry  of  thresholds.  Thresholds  increase  with  eccentricity  from 
the  fovea  starting  at  0.8°,  and  as  target  size  diminishes. 


Other  target  threshold  databases,  which  augment  the  VIDEM  model,  are  available 
in  the  model  when  the  visual  task  cannot  be  reduced  to  the  detection  of  disk-like  objects. 
Two  inside-the-cockpit  visual  tasks  were  included;  head-up  display  (HUD)  and  head- 
down  display  (HDD)  symbols.  The  latter  is  representative  of  emissive  CRT  displays. 
HUD  and  HDD  symbols  differ  slightly  in  their  computations  due  to  the  way  background 
light  is  treated.  The  contrast  ratio  is  formed  in  the  normal  way  (see  Equation  9)  for  HDD 
symbols,  but  the  background  light  is  added  back  into  the  numerator  for  HUD  symbols,  so 
the  contrast  ratio  is  simply: 


For  HUD  symbols  contrast  threshold  is  calculated  by; 

Ct  =  0.2667  •  e'(0-55S5*1°g(Bi))  (13) 

where  BL  =  Background  Luminance  in  cd-m'2.  For  HDD  symbols  the  contrast  threshold  is 
calculated  as: 


Ct  =  io(m*log(Si)+^ 


(14) 


where  m  and  b  take  on  values  given  in  Table  1  according  to  the  adapting  background 
luminance  level. 
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Table  1.  Slope  (m)  and  intercept  (b)  values  over  a  range  of  adapting  background 
luminances  for  determining  HDD  symbol  contrast  thresholds. 


f  (cd  m'2) 

§ 

<340 

<1700 
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-.44 
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-.01 

b 

-.29 

-.69 

-.93 

-1.16 

Figure  5  shows  how  HUD  and  HDD  contrast  thresholds  vary  with  adapting 
background  luminance  level.  At  low  levels  contrast  threshold  is  high,  but  declines  rapidly 
as  luminance  increases. 


Figure  5.  HUD  and  HDD  symbol  contrast  thresholds  as  a  function  of  adapting 
background  luminance  level.  Contrast  thresholds  are  high  in  the  dark,  but  decrease 
rapidly  as  background  luminance  level  increases.  Dawn  and  dusk  adapting  levels 
are  about  1  cd-m'2,  office  light  levels  are  about  100  cd-m'2,  and  bright  daylight  levels 
are  in  the  10000  cd-m'2  range. 

Wavelength  effects 

Because  lasers  are  monochromatic,  possible  wavelength  effects  on  glare  must  be 
considered.  Wavelengths  may  be  scattered  differently  by  the  ocular  media,  such  as  in 
Rayleigh  scattering.  Although  earlier  investigations  have  presented  contradictory  results, 
more  recent  research  indicates  that  entoptic  scatter  is  not  wavelength  dependent,  at  least 
when  extended  sources  are  used  (Wooten  &  Geri,  1987). 

Another  possible  effect  of  wavelength  on  glare  is  through  chromatic  contrast. 
Targets  of  equal  brightness  may  still  be  discriminable  by  a  difference  in  color  (Pokomy  & 
Smith,  1986).  Monochromatic  laser  light  might  interact  with  a  target’s  color  to  produce 
differential  effects  on  chromatic  sensitivity.  Indeed,  several  reports  have  shown 
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differential  effects  of  wavelength  on  colored  targets  (Previc,  1987;  Schmeisser,  1987; 
Varner  et  al.,  1988).  Responses  to  colored  targets  were  reduced  when  the  laser  and  the 
target  were  similar  in  color.  However,  these  effects  tended  to  be  small  compared  to 
impact  of  the  source  brightness  on  colored  target  visibility.  Consequently,  the  model  does 
not  incorporate  wavelengths  effects  except  for  correcting  the  source  irradianee  for  the 
chromatic  sensitivity  of  the  eye. 

LTMC  Flashblindness  Model 

Flashblindness  is  analyzed  the  same  way  as  glare,  except  that  the  Flashblindness 
Model  estimates  recovery  of  visual  sensitivity  after  the  source  has  been  turned  off.  As  a 
result,  the  magnitude  of  the  flashblindness  effect  depends  on  exposure  duration.  The 
Flashblindness  Model  analyzes  effects  for  single-pulse,  single  wavelength,  point  sources. 
However,  as  with  glare  sources,  multiple-pulse  exposures  greater  than  60  Hz  may  also  be 
analyzed  as  a  single  continuous  pulse. 

The  model  accepts  optical  radiation  inputs  in  terms  of  a  radiant  exposure,  He 
(J-cm'  )  and  exposure  duration  (s).  The  radiant  exposure  is  measured  at  the  cornea.  The 
radiant  exposure  is  converted  to  photometric  units  by  the  equation: 

Hv=6%5xHexVKxW,  (12) 

where  Hv  is  the  integrated  illuminance  in  hrm*2,s. 

The  Flashblindness  Model  adds  one  more  step  to  the  Glare  Model.  This  step 
represents  a  gradual  return  of  visual  sensitivity  after  the  optical  source  as  been  turned  off. 
Once  the  source  has  been  extinguished,  sensitivity  begins  to  recover  in  a  process  known  as 
dark  adaptation.  Just  as  the  Glare  Model,  the  Flashblindness  Model  uses  the  equivalent 
background  principle  to  model  the  recovery  in  sensitivity  (Crawford,  1947).  Flash  energy 
is  again  represented  as  a  veil  of  light,  but  now  one  that  fades  with  time.  The  initial  step  in 
the  Flashblindness  Model  is  the  expression  of  the  flash  energy  in  terms  of  the  integrated 
retinal  illuminance.  The  PSF  is  applied  to  the  flash  energy  to  estimate  the  retinal 
distribution  of  the  flash  as  an  integrated  luminance  in  cd-m'2-s,  that  is,  as  an  EBL  times  the 
exposure  duration,  t.  Then  the  integrated  luminance  is  expressed  in  terms  of  integrated 
retinal  illuminance,  Lv,  in  td-s.  This  is  done  by  multiplying  by  the  area  of  the  pupil,  and 
the  exposure  duration: 


Lv=PAxLegxt,  (13) 

where  Leq  is  the  EBL  in  cd-m'2,  t  is  the  exposure  duration  in  s,  and  PA  is  the  pupil  area  in 
mm2. 

Pupil  size 

Pupil  size  is  important  because  it  determines  the  amount  of  light  entering  the  eye, 
and,  therefore,  the  retinal  image  brightness.  It  is  also  important  because  flashblindness 
studies  quantify  flash  energy  in  terms  of  retinal  illumination,  so  that  input  into  the 
algorithm  must  be  expressed  in  td.  Pupil  size  is  normally  determined  by  the  ambient  level 
of  illumination.  The  flashblindness  source’s  duration  of  less  than  150  ms  would  be  too 
brief  to  elicit  a  sustained  pupillary  response.  Therefore,  the  adapting  background  drives 
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pupil  size  in  most  situations.  The  model  uses  an  algorithm  developed  by  de  Groot  and 
Gebhard  (1952)  to  estimate  pupil  size  from  adaptation  level.  The  relation  of  pupil 
diameter  to  the  brightness  of  an  adapting  light  source  is: 

Log  D  =  0.8558  -  0.0004(Log5  +  7.5972)3  (14) 

where  D  is  the  pupil  diameter  in  mm  and  B  is  the  adapting  field  brightness  in  cd-m'2. 

Figure  6  shows  the  association  between  pupil  diameter  and  ambient  luminance.  The 
equation  is  the  best  fit  to  eight  sets  of  data. 

The  model  limits  the  possible  range  of  pupil  diameters  between  2  and  7  mm. 
These  limits  were  chosen  because  the  adult  pupil  is  not  normally  less  than  2  mm  and  the 
ANSI  standard  for  the  estimation  of  laser  ocular  hazards  uses  a  maximum  of  7  mm  to 
compute  total  intraocular  energy  (American  National  Standards  Institute,  2000). 


Background  luminance  (log  cd  rri2) 


Figure  6.  Pupil  diameter  determined  by  ambient  adapting  luminance,  adapted  from 
de  Groot  and  Gebhard  (1952).  The  model  limits  pupil  diameter  between  2  and  7 
mm. 

Time-intensity  relationships  of  flash  brightness 

It  is  important  for  modeling  purposes  to  accurately  estimate  the  perceived 
brightness  of  the  flash.  The  flash  determines  the  initial  adaptation  level  and,  thus,  the 
initial  brightness  of  the  EBL.  As  with  glare,  this  can  be  done  by  comparing  the  flash 
illuminance  to  an  EBL  with  the  same  masking  effect.  However,  some  care  must  be  taken 
when  doing  this  because  the  brightness  of  a  flash  depends  on  intensity  and  duration  for 
brief  lights.  Common  observation  illustrates  that  the  brightness  of  a  steady  light  remains 
constant;  it  does  not  look  brighter  the  longer  one  stares  at  it.  However,  if  a  light  is  turned 
on  for  a  brief  time,  it  looks  dimmer  the  shorter  the  exposure  duration.  Therefore,  time 
must  be  considered  when  estimating  how  bright  a  light  flash  looks. 
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A  considerable  body  of  evidence  has  shown  that  the  visual  system  sums  the 
intensity  of  light  over  time  for  short  exposure  durations  (Hood  &  Finkelstein,  1986; 
Watson,  1986),  This  process  of  temporal  integration  is  known  as  Bloch’s  Law.  Bloch’s 
Law  states  that  the  perceived  brightness  of  a  light  is  the  product  of  its  intensity  and 
duration.  However,  this  law  is  valid  only  up  to  a  critical  length  of  time.  After  that,  no 
further  summation  takes  place  and  brightness  depends  on  intensity  alone.  This  critical 
length  of  time  is  known  as  the  critical  duration.  Bloch’s  Law  may  be  stated 
mathematically  as 


IxT  =  KforT<Tc,  (15) 

and 

I  =  KforT>Tc,  (16) 

where  I  =  intensity,  T=  time,  K  is  a  constant,  and  Tc  is  the  critical  duration.  This  law  is 
an  extension  of  the  Bunsen-Roscoe  law  of  photochemistry.  A  schematic  diagram  of 
Bloch’s  law  is  illustrated  in  Figure  7.  In  log-log  coordinates,  Equation  (15)  gives  a 
straight  line  with  a  slope  of  -1. 


Figure  7.  An  idealized  threshold-duration  function  (Bloch’s  law).  The  function 
describes  thresholds  for  rectangular  light  pulses  as  a  function  of  duration.  The  left 
limb  of  the  function  has  a  slope  of — 1  and  the  right  limb  has  a  slope  of  0.  The 
transition  between  the  two  limbs  occurs  at  the  critical  duration  Tcand  the  critical 
intensity  Ic.  The  inset  shows  the  wave  form  of  a  rectangular  pulse  with  duration  T 
and  intensity  I.  From  Watson  (1986). 

One  consequence  of  this  law  is  that  lights  of  equal  energy  are  equally  detectable, 
that  is,  intensity  and  duration  are  reciprocal.  Any  combination  of  intensity  and  duration 
that  yields  a  constant  product  appears  equal  in  brightness.  The  intensity  must  be  above  the 
critical  intensity  and  the  duration  must  be  below  the  critical  duration.  Equation  (16)  yields 
a  horizontal  line  in  Figure  7  and  indicates  that,  after  the  critical  duration,  brightness 
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depends  on  intensity  alone.  When  experimental  data  are  considered,  rarely  do  the  two 
functions  show  such  an  abrupt  shift  from  one  slope  to  the  other.  Often  a  transition  region 
is  observed  from  complete  reciprocity  to  no  reciprocity.  This  region  of  partial  summation 
is  sometimes  reported  as  having  a -0.5  slope,  and  so  follows  an  inverse-square  law. 

A  number  of  studies  have  tried  to  pin  down  the  exact  length  of  the  critical  duration 
and  have  found  that  it  is  influenced  by  many  factors  including  light  adaptation  level 
(Barlow,  1958;  Crawford,  1937;  Stewart,  1972),  the  size  of  the  light  (Barlow*  1958;  Hood 
&  Finkelstein,  1986;  Kam,  1936),  and  wavelength  (Sperling  &  Jolliffe,  1965).  The 
shortest  period  of  time  for  which  reciprocity  holds  was  investigated  by  Brindley  (1952). 
He  investigated  equal  energy  suprathreshold  lights  between  0.4  /is  and  8.9  ms  and  found 
that  reciprocity  was  valid  down  to  the  shortest  duration.  Brindley  was  limited  in  the  total 
energy  he  could  deliver  at  the  shortest  durations  because  of  the  limited  intensity  of  his 
light  sources.  As  a  result  the  0.4  /is  light  was  not  very  bright,  only  0.3  td  s.  It  would  be 
beneficial  to  re-examine  the  reciprocity  issue  with  brighter  sources.  Today’s  lasers  with 
high  peak  power  and  short  pulsewidths  would  be  ideal  to  investigate  the  temporal 
integration  properties  of  the  visual  system  at  heretofore-untested  durations  and  intensities. 

As  a  general  rule  one  may  conclude  that  reciprocity  is  valid  up  to  a  maximum  time 
of  about  150  ms  in  the  cone  system.  Partial  reciprocity  occurs  from  150-1000  ms  and 
brightness  is  largely  independent  of  duration  after  1  s  (Kosnik,  1998;  Watson,  1986).  The 
following  rule  has  been  adopted  for  flashblindness  modeling.  The  brightness  of  a  flash  is 
proportional  to  its  energy  for  durations  of  150  ms  and  less. 

Time-intensity  relationships  in  flashblindness  recovery 

Although  reciprocity  for  brightness  has  been  fairly  well  established,  reciprocity  in 
flashblindness  recovery  is  not  nearly  so  well  understood.  Consequently,  a  valid 
flashblindness  model  must  establish  some  limits  of  reciprocity  where  equivalent  energy 
exposures  have  the  same  effect  on  dark  adaptation.  Like  brightness  perception,  the  bulk 
of  the  evidence  indicates  reciprocity  for  short-duration  exposures  between  0.1  and  8.5  ms 
(Bowie  &  Collyer,  1973;  Chisum,  1973;  Miller,  1965).  However,  studies  comparing 
millisecond  and  second  durations  generally  did  not  find  reciprocity  (Bowie  &  Collyer, 
1973;  Crawford,  1946;  Miller,  King,  &  Schloessler,  1968).  The  evidence  for  durations 
between  10  ms  and  1  s  is  mixed.  Crawford  (1946)  found  that  durations  between  9  and  900 
ms  showed  reciprocity,  more  so  for  lower  energy  exposures  than  for  higher  ones.  Bowie 
and  Collyer  (1973)  found  only  partial  reciprocity  for  durations  between  13  and  100  ms 
and  no  reciprocity  between  100  ms  and  1  s.  Thus,  it  appears  that  reciprocity  is  only  found 
under  a  restricted  set  of  conditions.  For  purposes  of  the  model,  therefore,  equal-energy 
exposures  having  durations  less  than  or  equal  tol50  ms  are  modeled  with  the  same  dark 
adaptation  function.  Although  the  literature  is  not  in  complete  agreement  on  the  higher  end 
of  this  range,  150  ms  is  a  reasonable  limit  because  the  eye  blink  reflex  at  150  ms  would 
prevent  a  continuous  duration  exposure  from  occurring  anyway  (Irving,  1960). 
Furthermore,  eye  movements  would  come  into  play  after  150  ms,  further  distorting  the 
shape  of  the  retinal  light  distribution  (Kosnik,  1988;  Stamper,  Lund,  Molchaney,  &  Stuck, 
2000).  Therefore,  the  model  handles  most  combinations  of  duration  and  intensity  that 
would  likely  occur.  Nevertheless,  further  research  should  be  conducted  to  firmly  establish 
the  applicability  of  the  reciprocity  law  in  the  range  of  10  ms  to  1  s. 
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The  reciprocity  rule  is  extended  to  include  the  integration  of  multiple  pulses.  Any 
train  of  pulses  within  150  ms  may  be  considered  a  single  continuous  pulse  as  far  as 
flashblindness  modeling  is  concerned.  The  individual  pulses  of  a  multiple  pulse  train  are 
simply  summed  together.  The  flashblindness  effect  is  considered  the  same  no  matter  how 
the  pulses  are  distributed.  Furthermore,  the  model  assumes  that  the  visual  system  has 
reached  a  steady-state  level  of  adaptation.  Thus,  it  does  not  account  for  changing 
adaptation  effects  during  the  initial  pulses  of  a  multiple  pulse  train. 

The  Flashblindness  Model  also  has  an  upper  exposure  limit.  Based  on  the 
experimental  evidence,  which  suggests  that  if  all  of  the  photopigment  is  bleached  by  a 
given  exposure,  then  any  additional  energy  has  no  effect  on  the  light  adapted  state.  The 
total  energy  of  7.6  log  td  s  was  set  as  the  upper  limit  as  this  exposure  bleaches  all  of  the 
photopigment  (Rushton  &  Henry,  1968).  This  limit  does  not  imply  that  higher  energy 
flashes  have  no  effect  on  other  aspects  of  the  model.  For  example,  the  retinal  light 
distribution  changes  with  higher  energy  flashes,  and  this  aspect  is  modeled  in  the  PSF. 


Flashblindness  Recovery  Functions 

AFRL/HEDO  interest  in  visual  recovery  centers  around  exposure  to  lasers  under 
operational  conditions.  Such  exposures  are  de  facto  unintentional  and  it  is  safe  to  assume 
that  the  individual  will  make  every  effort  to  limit  the  exposure  under  most  circumstances. 
Therefore,  the  model  focuses  on  brief  exposures.  Visual  recovery  immediately  following 
laser  exposure  is  a  primary  concern  to  AFRL/HEDO.  Air  Force  personnel  most  at  risk  for 
laser  exposure  are  those  engaged  in  critical  tasks,  such  as  flying  an  aircraft.  The  ability  to 
maintain  continuous  flight  control  is  of  utmost  importance.  For  this  reason,  visual 
recovery  within  the  first  two  minutes  is  the  focus  of  the  modeling  effort. 

Three  flashblindness  recovery  functions  were  considered  for  use  in  LTMC.  The 
functions  all  have  the  same  general  form  in  that  they  express  flashblindness  recovery  in 
terms  of  an  equivalent  background  luminance  (EBL),  which  decays  as  a  function  of  time 
after  the  flash.  Previous  AFRL/HEDO  flashblindness  models  have  used  a  flashblindness 
recovery  function  developed  by  Menendez  and  Smith  (1990).  Smith  later  carried  out 
extensive  studies  of  flashblindness  recovery  in  human  subjects  and  developed  a  modified 
function  (Smith,  1996).  In  addition,  an  exponential  decay  function  is  considered  (Hahn  & 
Geisler,  1995). 

Just  as  in  glare  modeling,  the  equivalent  background  principle  is  a  useful  construct 
for  describing  the  light  adapted  state  of  the  eye;  for  flashblindness,  the  EBL  includes  a 
dark  adaptation  component  as  well  (Crawford,  1947).  Some  have  argued  that  the  EBL 
represents  the  true  state  of  adaptation  at  any  time  during  dark  adaptation  (Crawford,  1947; 
Miller,  1965, 1966a,  1966b),  but,  this  idea  is  still  controversial  (Hahn  &  Geisler,  1995; 
Hood  &  Finkelstein,  1986;  Kosnik,  1998).6  One  advantage  of  using  the  EBL  for  modeling 


“The  EBL  has  been  compared  to  a  flash  afterimage,  which  has  also  been  used  to  represent  the  light 
adapted  state  of  the  eye  (Miller,  1966a,  1966b).  The  EBL  can  act  like  an  afterimage  under  special 
circumstances.  Both  decrease  in  brightness  with  post-exposure  time.  However,  the  two  are  not  the  same 
(Barlow  &  Sparrock,  1964).  The  EBL,  on  the  one  hand,  is  a  hypothetical  construct.  An  afterimage,  on  the 
other  hand,  is  a  physiological  phenomenon  that  has  a  perceptual  counterpart.  It  can  be  positive  or  negative 
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flashblindness  recovery  is  that  it  is  independent  of  target  and  background  conditions,  so 
the  flash  effect  can  be  generalized  across  all  viewing  conditions. 

The  Menendez-Smith  (1990)  function  and  the  Smith  (1996)  function  have  the 
same  general  form.  The  functions  relate  the  EBL  magnitude  to  the  initial  flash  energy  as  a 
function  of  time  after  the  flash.  Each  function  is  a  log-log  equation  of  the  form: 


Log  EBL  =  m  Log  t  +  c,  (17) 

where  EBL  is  in  cd  m'2  and  t  is  time  after  the  flash  in  s.  The  slope  of  the  function,  m, 
determines  the  rate  at  which  the  EBL  fades  and  c  is  the  offset.  The  difference  between 
the  two  models  is  how  they  are  parameterized. 

Hie  third  function  is  based  on  the  Hahn-Geisler  (Hahn  &  Geisler,  1995)  dark 
adaptation  model.  This  model  uses  an  exponential  decay  function  to  describe  the  dark 
adaptation  process.  It  is  based  on  the  work  of  Rushton  (1964;  Rushton  &  Henry,  1968) 
and  Hollins  and  Alpem  (1973).  The  model  was  adapted  to  make  it  into  the  form  of  an 
EBL.  Its  general  form  is: 

Log(£5L  +  l)  =  ouf'"',  08) 


where  a  is  the  initial  light  adaptation  level  at  flash  offset,  t  is  time  in  s,  and  tc  is  the  time 
constant  of  decay.  The  EBL  is  continuous  from  the  offset  of  the  flash  and  decays  to  dark 
in  this  function.  The  forms  of  the  three  functions  are  shown  in  Figure  8  for  a  flash  of  7.6 
log  td  s.  Note  that  all  three  functions  convert  retinal  flash  energy  into  an  external  EBL. 


Figure  8.  The  forms  of  the  three  flashblindness  recovery  functions  for  an  input  flash 
energy  of  7.6  log  td*s. 


(white  or  gray),  or  colored,  and  it  can  change  in  polarity  and  color  with  time.  Its  appearance  depends  on  the 
background  on  which  it  is  seen  and  it  can  be  refreshed  with  eye  blinks  or  eye  movements. 
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Comparison  of  the  flashblindness  functions 

The  three  functions  were  compared  against  a  set  of  flashblindness  data  to 
determine  which  function  best  represented  the  flashblindness  recovery  process.  First, 
though,  the  functions  had  to  be  parameterized.  The  Menendez-Smith  function  has  two 
free  parameters,  m  and  c.  The  slope  of  the  decay  function,  m,  was  assigned  a  value  of -3. 
This  quantity  is  based  on  two  studies  by  Miller  (Miller,  1965, 1966b)  who  found  that  the 
log  brightness  of  an  afterimage  faded  at  a  rate  of  -3  log  t.  Miller  found  that  this  rate  was 
constant  for  flash  energies  between  6.68  and  7.48  log  td-s.  Menendez  and  Smith  (1990) 
generalized  this  rate  to  all  flash  energies. 

The  offset,  c,  was  determined  by  noting  that  c  =  3Log  t  +  0.5  when  the  EBL 
decayed  to  3.2  cd-m'2.  It  was  assumed  that  the  target,  a  Sloan-Snellen  letter,  0.22  cd-m'2  in 
brightness,  was  just  resolvable  when  the  EBL  reached  this  level.  In  other  words,  the  letter 
had  a  putative  resolution  threshold  of  7%.  Menendez  and  Smith  (1990)  used  the  results  of 
Miller’s  (1965)  flashblindness  recovery  data  to  determine  at  what  time  the  EBL  faded  to 
the  extent  that  the  target  reached  threshold.  Miller  reported  in  her  1965  study  five 
recovery  times  based  on  five  different  flash  energies.  By  substituting  the  recovery  times 
into  Equation  (17),  Menendez  and  Smith  (1990)  derived  five  values  for  c.  These  values 
were  plotted  against  flash  energy  to  find  a  general  expression  for  c.  The  results  are  shown 
in  Figure  9.  It  can  be  seen  that  the  c  values  form  a  straight  line  when  plotted  on  a  log-log 
chart.  The  best  fitting  line,  fit  by  the  least  squares  method,  was: 

c  =  1.75F-6.33,  (19) 

where  F  is  the  flash  energy  in  log  td-s. 


Figure  9.  Derivation  of  the  offset  parameter,  c,  in  the  Menendez-Smith  function. 

The  coefficient,  c,  was  derived  from  five  recovery  trials  collected  at  five  different 
flash  energies.  A  general  equation  for  c  was  obtained  by  fitting  a  linear  function. 

The  coefficient  of  determination,  R2,  was  0.98,  which  indicated  that  the  nearly  all  of  the 
variation  in  the  data  was  accounted  for  by  the  regression.  Substituting  for  c,  the  decay 
function  became: 
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(20) 


Log(EBL)=  -3Log(f)+  1.75F  -  6.33 , 

The  function  was  modified  to  limit  flash  energy  to  7.6  log  tds  and,  thus,  became: 

Log(EBL)  =  -3Log(f)+ 1.75  x  min(7.6,  F)  -  6.33 .  (21) 

The  Smith  and  Hahn-Geisler  functions  were  parameterized  using  laser 
flashblindness  recovery  data  from  Smith  (1996).  This  data  set  was  chosen  because  it 
collected  recovery  times  over  a  wide  range  of  flash  energies.  Thus,  the  form  of  visual 
recovery  could  be  assessed  over  a  wide  dynamic  range.  Recovery  times  were  measured  in 
two  observers  for  identification  of  a  90%  contrast  Landolt  C  at  seven  background 
luminance  levels  ranging  from  0.34  to  154  cd  m'2.  Ten  laser  flash  energies  were  used 
ranging  from  5.41  to  7.57  log  td-s.  The  laser  exposures  were  generated  from  the  514-nm 
line  of  a  Coherent  argon  ion  laser  and  were  100-ms  in  duration.  Each  exposure  was 
presented  as  a  circular,  uniform  Maxwellian-view  spot  subtending  1.15°  at  the  eye.  In 
addition,  veiling  luminance  thresholds  were  measured  at  each  of  the  background 
luminance  levels  and  used  to  transform  the  recovery  data  into  EBL  versus  time.  These 
data,  shown  in  Appendix  A,  were  used  to  parameterize  the  functions. 

The  two  free  parameters  in  Equation  (17),  m  and  c,  were  fit  simultaneously  using 
the  least  squares  method  of  linear  regression.  The  solver  function  of  Microsoft  Excel™ 
performed  the  fitting  procedure.  Solver  successfully  converged  on  a  solution  in  each  case. 
The  resulting  estimated  parameter  values  are  shown  in  Table  2  for  each  flash  energy. 

The  same  procedure  was  used  to  parameterize  the  Hahn-Geisler  function.  Solver 
was  used  again  to  simultaneously  fit  the  two  free  parameters,  a  and  tc,  in  Equation  (18). 
Solver  successfully  converged  on  a  solution  and  the  resulting  estimated  values  are  shown 
in  Table  2. 

Table  2.  Empirically  derived  free  parameter  estimates  of  the  Smith  and  Hahn- 

Geisler  functions. 
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Both  functions  fit  the  data  well.  All  the  residual  errors  from  the  least  squares 
regression  procedures  were  small.  An  example  of  the  fits  of  the  two  functions  for  the  7,57 
log  td*s  flash  is  shown  in  Figure  10.  The  results  of  the  fitting  procedure  for  the  other  flash 
energies  are  shown  in  Appendix  B. 

The  fitting  exercise  showed  that  both  functions  simulated  the  flashblindness 
recovery  process  well.  Nevertheless,  it  would  be  advantageous  to  find  a  general  function 
that  does  not  depend  on  data  to  estimate  the  free  parameters.  Thus,  an  attempt  was  made 
to  express  the  free  parameters  in  terms  of  flash  energy.  This  was  done  by  plotting  die  free 


Figure  10.  Results  of  the  parameter  estimation  exercise  using  the  Smith  (1996)  data. 

The  plot  shows  the  data  from  the ,7.57  log  td*s  flash.  The  data  are  shown  as  the 
diamond  symbols;  the  solid  line  is  the  Smith  function  and  the  broken  line  is  the 
Hahn-Geisler  function. 

parameter  values  obtained  in  the  data  fitting  exercise  against  flash  energy  and  examining 
the  results  for  trends.  A  least-squares  regression  procedure  was  used  to  derive  equations 
for  Smith’s  parameters,  m  and  c.  The  parameter  values  obtained  from  the  data  fitting 
procedure  are  shown  in  Figure  1 1 . 

The  squares  indicate  the  slope  parameter,  m  and  the  diamonds  and  circles  show  the 
offset  parameter,  c.  The  m  term  was  fit  by  two  straight  lines.  (See  Smith,  1996,  for  the 
rationale  for  this  approach.)  The  first  four  data  points  were  fit  by  a  straight  line  (R2  = 
0.71)  of  the  form: 


m  =  F-8.92,  F  <  6.11  logtds  (22) 

The  rest  of  the  data  points  were  fit  by  the  line  (R2  =  0.03): 

m  =  -0.022F  -  2.73,  F>6.11  logtds.  (23) 

However,  the  slope  of  this  line  was  not  significantly  different  from  zero  so  the  equation 
reduced  to  the  constant  term,  -2.73.  The  second  parameter,  c,  was  well  fit  (R2  =  0.99)  by 
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(24) 


a  second-order  polynomial  of  the  form: 

c  =  6.4F  -  0.41F2  - 19.43 . 
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Figure  11.  Parameter  estimates  obtained  from  fitting  the  Smith  function  to  Smith’s 
(1996)  flashblindness  data.  The  parameter  estimates  are  plotted  against  flash 
energy.  The  offset  parameter,  c,  is  shown  as  the  square  symbols  and  uses  the  left 
ordinate.  The  diamonds  and  circles  indicate  the  slope  parameter,  m,  which  uses  the 
right  ordinate.  The  three  lines  show  the  results  of  the  curve-fitting  procedures.  The 
two  solid  lines  show  the  linear  trends  for  the  m  parameter  and  the  dashed  line  shows 
the  second  order  polynomial  fit  of  the  c  parameter. 


The  plots  of  the  two  Hahn-Geisler  parameters,  a  and  tc,  are  shown  in  Figure  12. 
The  diamonds  and  circles  show  the  a  parameter  and  the  triangles  indicate  the  time 
constant  parameter,  tc.  As  shown  by  Smith  (1996),  the  a  term  was  fit  by  two  straight  lines. 
The  first  three  data  points  were  fit  by  a  straight  line  (R2  =  1.0)  of  the  form: 


a  =  -1.49F  +  12.19,  F  <  6.0  log  td-s .  (25) 

The  rest  of  the  data  points  were  fit  with  a  line  (R2  =  0.54)  of  the  form: 

a  -  0.106F  +  2.68,  F  £  6.0  log  td-s .  (26) 


The  second  parameter,  tc,  was  fit  by  a  cumulative  Gaussian  function  of  the  form: 


(27) 

s 


where  k  is  a  constant  and  p  and  c  are  the  mean  and  standard  deviation  of  the  distribution. 
This  function  was  fit  to  the  data  by  Solver.  The  function  fit  the  data  best  (R2  =  0.97)  when 
k  =  43.23,  p  =  6.23,  and  a  =  0.60.  The  cumulative  Gaussian  function  seemed  an 
appropriate  model  because  it  fit  the  limits  of  the  flash  energy.  Flash  energy  cannot  go 
below  zero  and  it  asymptotes  after  the  photoreceptors  become  saturated,  at  about  7.6  log 
td-s  (Rushton  &  Henry,  1968). 
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Figure  12,  Parameter  estimates  obtained  from  fitting  the  Hahn-Geisler  function  to 
Smithes  (1996)  fiashbiindness  data.  The  parameter  estimates  are  plotted  against 
flash  energy.  The  a  parameter,  scaled  by  the  left  ordinate,  is  shown  as  the 
diamonds  and  circles.  The  triangles  represent  the  time  constant  parameter,  tc, 
which  is  scaled  by  the  right  ordinate.  The  three  lines  show'  the  results  of  the  curve 
fitting  procedures.  The  two  solid  lines  show  the  linear  trends  for  the  a  parameter 
and  the  dashed  line  shows  the  cumulative  Gaussian  fit  of  the  time  constant. 


The  flash  energies  were  substituted  back  into  the  Equations  (22-24)  for  the  Smith  function 
and  Equations  (25-27)  for  the  Hahn-Geisler  function  to  obtain  model  parameters  at  each 
flash  energy.  Table  3  shows  the  predicted  model  parameters. 

Table  3.  Predicted  estimates  of  the  free  parameters  in  the  Smith  and  Hahn-Geisler 
functions. 
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A  Test  of  the  Flashblindness  Recovery  Functions 

The  three  functions,  with  parameters  derived  from  flash  energy,  were  fit  again  to 
the  Smith  (1996)  data  and  one  other  data  set  from  a  study  by  Menendez  and  Garcia 
(1985).  The  latter  study  was  chosen  because  it  employed  a  similar  methodology  to  the 
Smith  (1996)  study.  Menendez  and  Garcia  used  a  single  flash  energy  of  6.7  log  tds,  an 
energy  also  used  by  Smith  (1996),  and  both  used  a  uniform,  spatially  extended  source,  and 
a  Maxwellian-view  system.  It  would  be  expected  that  the  decay  functions  from  both 
studies  would  have  the  same  time  course. 

The  Menendez-Smith,  Smith,  and  Hahn-Geisler  functions  were  plotted  against  the 
Smith  (1996)  flashblindness  recovery  data.  Some  of  these  plots  are  shown  in  Figure  13. 
Each  plot  represents  one  flash  energy.  The  fits  to  the  other  flash  energies  are  shown  in 
Appendix  C.  The  Smith  and  the  Hahn-Geisler  functions  provided  good  fits  to  the  Smith 
(1996)  data.  The  Menendez-Smith  function,  on  the  other  hand,  frequently  overestimated 
or  underestimated  the  data.  This  was  probably  due  to  the  use  of  a  log-linear  function  to 
estimate  the  offset  parameter,  c,  which  means  that  it  continued  to  increase  with  flash 
energy  (see  Figure  9).  The  parameters  of  the  other  two  functions  asymptote  at  higher 
flash  energies  (see  Figure  1 1  and  Figure  12).  The  asymptotic  behavior  of  these  parameters 
follows  from  the  limit  of  flash  energy  on  photopigment  bleaching.  The  parameters  should 
start  to  asymptote  as  flash  energy  approaches  7.6  log  td-s  as  no  further  bleaching  occurs 
and  the  light-adapted  state  of  the  retina  has  reached  a  maximum. 

The  three  functions  were  plotted  against  data  from  the  Menendez  and  Garcia 
(1985)  study.  These  data,  shown  in  Figure  14,  allowed  for  an  unbiased  test  of  the 
functions  because  they  were  not  used  in  parameterization  processes.  The  Smith  and 
Hahn-Geisler  functions  fit  the  best.  The  Menendez-Smith  function  tended  to 
underestimate  recovery  times,  particularly  when  the  EBL  decayed  to  a  very  low  level. 

Based  on  the  foregoing  tests,  the  performances  of  the  Smith  and  Hahn-Geisler 
functions  were  too  similar  to  pick  one  function  over  the  other.  However,  upon  further 
examination  it  became  apparent  that  one  function  performed  better  than  the  other.  The 
Smith  function  was  more  robust  in  that  it  readily  extrapolated  to  flash  energies  outside  the 
range  used  to  parameterize  the  function.  For  example,  Smith  (1996)  investigated 
flashblindness  recovery  to  lower  flash  energies  in  a  chopped  pulse  experiment.  He  found 
that  the  Smith  function  accurately  predicted  up  to  1-s  of  flashblindness,  however,  when  we 
tested  the  Hahn-Geisler  model,  it  did  not  predict  any. 

Furthermore,  the  Hahn-Geisler  function  had  some  trouble  predicting  recovery 
times  beyond  the  flash  energy  range  for  which  the  function  was  parameterized.  The 
reason  for  this  lack  of  generality  can  be  found  in  the  equation  used  to  parameterize  the 
time  constant  (see  Figure  12).  The  time  constant,  tc,  tails  off  to  zero  by  the  time  the  flash 
energy  reaches  4.5  log  td  s.  It  is  clear  that  the  tc  tapers  off  too  quickly,  as  it  should  not 
have  reached  zero  until  the  exposure  reached  0  log  td-s.  This  behavior  of  the  Gaussian 
function  was  probably  to  due  the  lack  of  data  points  at  lower  flash  energies.  If  some  data 
were  collected  at  flash  energies  below  5.4  log  td  s,  then  the  cumulative  Gaussian  would 
have  tapered  off  more  gradually.  When  future  studies  provide  data  at  lower  flash  energies, 
the  Hahn-Geisler  function  may  yet  prove  to  be  a  viable  predictor  of  flashblindness 
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Figure  13,  The  fit  of  the  three  functions  to  the  Smith  (1996)  flashblindness  recovery 
data  for  three  flash  energies.  Flash  energy  is  noted  at  the  top  of  each  plot.  The 
t  Smith  is  the  solid  line,  the  Hahn-Geisler  function  is  the  broken  line  and  the 

Menendez-Smith  function  is  the  dotted  line. 
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recovery.  At  the  present  time,  the  Smith  function  provided  the  most  accurate  predictions 
across  a  wide  range  of  flash  energies  and  is,  therefore,  adopted  for  LTMC. 
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Figure  14.  The  three  decay  functions  are  plotted  against  the  Menendez  and  Garcia 
(1985)  flashblindness  recovery  data.  The  Smith  (solid  line)  and  Hahn-Geisler 
(broken  line)  functions  give  the  best  fits.  The  Menendez-Smith  (dotted  line)  function 
underestimates  recovery  time  compared  to  the  other  two,  especially  at  low  EBLs. 


Summary 

This  report  described  the  status  of  flashblindness  and  glare  modeling  at 
AFRL/HEDO.  The  rationale  guiding  the  models’  development  and  the  underlying 
mathematical  functions  were  presented  and  were  supported  by  the  pertinent  scientific 
literature.  These  models  predict  the  size  and/or  duration  of  temporary  visual  scotomas 
resulting  from  optical  radiation  exposure.  The  basic  components  of  die  models  were 
described.  Those  components  included  the  wavelength,  pulse-width,  and  irradiance  of  the 
laser  source,  the  target  and  background  parameters,  and  the  optical  and  adaptation  state  of 
the  eye.  Three  candidate  flashblindness  recovery  functions  were  reviewed.  The  candidate 
functions  were  tested  against  an  independent  flashblindness  recovery  data  set  to  determine 
which  function  best  fit  the  available  data.  A  function  by  Smith  (1996)  was  found  to  be  the 
best  performing  function. 

The  following  specifications  are  adopted  for  the  LTMC  Glare  and  Flashblindness 
Models:  » 

•  Glare  is  defined  as  an  exposure  that  appears  continuous  for  at  least  one  second. 

•  The  maximum  light  adaptation* effect  of  the  glare  source  is  limited  to  7.6  log  td.  * 

•  Exposure  intensity  alone  determines  perceived  brightness  for  exposure  durations 
beyond  one  second. 
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•  The  Flashblindness  Model  is  limited  to  analyzing  exposures  150  ms  or  less. 

•  The  maximum  light  adaptation  effect  of  a  flashblinding  source  is  limited  to  7.6  log 
td-s. 

•  The  Flashblindness  Model  is  limited  to  analyzing  single-pulse  exposures  (or 
multiple  pulses  within  150  ms). 

•  The  Flashblindfiess  Model  uses  the  Smith  (1996)  function  for  predicting  the  time 
course  of  flashblindness  recovery. 

•  The  minimum  predicted  flashblindness  recovery  time  is  one  second. 

•  Exposure  energy  and  EBL  brightness  follow  the  reciprocity  law  for  exposures  up 
to  150  ms. 

•  Both  Models  are  limited  to  analyzing  point  source  exposures. 

•  Both  Models  are  limited  to  analyzing  single  wavelength  exposures. 

•  Both  models  use  the  TNO  point  spread  function  to  estimate  intraocular  light 
scatter. 

•  Both  models  use  the  photopic  spectral  luminosity  function  to  correct  for  the 
wavelength  sensitivity  of  the  eye. 

Recommendations 

The  following  recommendations  are  provided  to  guide  future  modeling  efforts: 

•  The  TNO  algorithm  for  estimating  the  retinal  light  distribution  of  extended  light 
sources  should  be  included  to  account  for  light  scattering  by  optical  media  such  as 
canopies  and  windscreens. 

•  The  Flashblindness  Model  is  designed  to  analyze  point  sources,  but  the 
flashblindness  recovery  function  was  developed  and  validated  on  data  from 
extended  source  exposures.  The  flashblindness  recovery  function  should  be  re¬ 
validated  using  point  source  exposures. 

•  A  multiple  pulse  model  should  be  developed  to  predict  the  effect  of  cumulative 
pulses  on  the  light  adaptation  level  of  the  visual  system. 

•  Experiments  should  be  conducted  to  resolve  reciprocity  issues  for  exposures 
between  10  ms  and  one  second. 

•  The  model  should  be  enhanced  to  allow  independent  positioning  of  the  target  and 
light  source  to  fall  anywhere  on  the  retina. 

•  The  model  should  be  enhanced  by  adding  targets  to  the  database  that  involve 
critical  flying  tasks,  especially  those  involving  HUD  symbology. 
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Appendix  A 

Recovery  time  (s)  data  versus  the  EBL  (cd  m2)  at  ten  laser  energy  levels  (log  td  s) 
from  Smith  (1996). 
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Appendix  B 

Plots  of  the  Smith  (solid  line)  and  Hahn-Geisler  (broken  line)  functions  fitted  to 
the  Smith  (1996)  data  (symbols;  see  Appendix  A),  using  the  empirically  derived 
parameters.  The  data  are  given  in  recovery  time  versus  the  EBL.  Flash  energies  in  log 
td*s  are  reported  at  the  top  of  each  figure. 
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Appendix  C 

Plots  of  the  Smith  (solid  line),  Hahn-Geisler  (broken  line),  and  Menendez-Smith 
(dotted  line)  functions  fitted  to  the  Smith,  1996,  data  using  the  parameters  derived  from 
the  flash  energies.  The  data  are  given  in  recovery  time  versus  die  EBL.  Flash  energies  in 
log  td-s  are  reported  at  the  top  of  each  figure. 
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Appendix  D 

Data  from  Menendez  and  Garcia  (1985)  showing  the  recovery  time  vs.  the  EBL. 
The  raw  data  were  transformed  into  EBLs  from  contrast  scores.  Data  are  from  four 
observers  and  five  sine  wave  gratings,  1,4,8,12,  and  24  c-d'1.  Each  score  represents  the 
average  from  ten  flash  recovery  trials. 
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